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We explore the cosmological solutions of classes of non-linear bigravity theories. These 
theories are defined by effective four- dimensional Lagrangians describing the coupled dy- 
namics of two metric tensors, and containing, in the linearized limit, both a massless gravi- 
ton and an ultralight one. We focus on two paradigmatic cases: the case where the coupling 
between the two metrics is given by a Pauli-Fierz-type mass potential, and the case where 
this coupling derives from five-dimensional brane constructions. We find that cosmological 
evolutions in bigravity theories can be described in terms of the dynamics of two "relativis- 
tic particles" , moving in a curved Lorenzian space, and connected by some type of nonlinear 
"spring" . Classes of bigravity cosmological evolutions exhibit a "locking" mechanism under 
which the two metrics ultimately stabilize in a bi-de-Sitter configuration, with relative (con- 
stant) expansion rates. In the absence of matter, we find that a generic feature of bigravity 
cosmologies is to exhibit a period of cosmic acceleration. This leads us to propose bigrav- 
ity as a source of a new type of dark energy ("tensor quintessence"), exhibiting specific 
anisotropic features. Bigravity could also have been the source of primordial inflation. 
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1 Introduction 



Recently, there has been some interest in multigravity theories |T],|2|,|3|,H,|5|,|6| where gravity is 
modified at cosmological scales. These theories involve brane configurations in a higher than 
four dimensional spacetime where normal four dimensional gravity on the branes is modified 
in the far infrared due to the presence of a massive (but ultralight) graviton component 
in the low energy theory (see @ for a review and || for a detailed presentation). The 
attractive feature of these theories is that they provide an alternative observational window 



to extra dimensional physics which is testable in current observations 0. [10]. [□], [12], |H| 



Most importantly, these modifications of gravity are at such scales from which current 



observations indicate a dark energy component in our universe []14||L5||. It would be tempting 
to attribute this to the dynamics of a multigravity system. 

The basic idea in constructing multigravity models is to localize gravity at the same 
time in different places along the extra dimension(s). Once one has a higher dimensional 
brane configuration which localizes gravity, e.g. Ul6|,|T7|l, the low energy effective theory is 



governed by a massless graviton field. By multilocalizing gravity |18| in a superposition 
of such configurations, the degeneracy of the massless modes will be lifted and the low 
energy theory will contain apart of a massless mode, a collection of light massive gravitons. 
How light these gravitons will be depends on how strong the localization in the single 
graviton configuration is. In the particular models that have been examined in the literature 
|j],0,[|,[|,|5|,|6|], the localization was exponentially strong and thus the mass splittings between 
the light gravitons and the remaining of the Kaluza-Klein (KK) spectrum was exponentially 
large. This gave the opportunity to realize models which escaped observational bounds 
and had interesting phenomenological implications. An alternative mechanism which has 
similar effects arises in the case where four dimensional gravity is induced on the brane 
due to quantum loops of matter living on the brane |l9|,p0|1 . In that case, the old result of 
Sakharov [21]] (for a review of induced gravity see [p2]) was exploited to render the geometry 
on a brane, embedded in a flat higher dimensional space, four dimensional. 

However, so far the models Jl], 0, [3], |j, [5], || where studied at the linearized level which 
becomes invalid when speaking about cosmological distance dynamics. One clearly has 
to go beyond the linear theory, which is the main aim of this paper. Nonlinear bigravity 
theories were first introduced in the seventies as effective descriptions of a sector of hadronic 
physics [p3 |. It is argued in a companion paper [24] that nonlinear bigravity theories can 
arise in several different (purely gravitational) contexts: multibrane configurations, certain 
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classes of Kaluza-Klein models, some types of non-commutative geometry models, etc. It 
is, therefore, important to try to delineate what are the generic predictions of classes of 
bigravity theories. One of the main conclusions of the present paper is that bigravity 
naturally gives rise to a late period of cosmic acceleration. Bigravity can then be used as 
a new theoretical model of dark energy (with specific anisotropic features, in certain cases, 
that make it phenomenologically distinguishable from quintessence models). 

Accelerating solutions were also found in the context of one particular model of brane- 



induced gravity [19,20], and their phenomenological consequences have been explored in 



detail as a possible theoretical model of dark energy P3, 2^|. Both models share the common 
feature of modifications of gravity at large scales. However, our work differs in several ways. 
We study general classes of four- dimensional effective theories instead of one particular five- 
dimensional model. We study general classes of solutions of these theories including their 
stability. We find several different types of accelerating solutions some being isotropic, but 
many featuring a novel type of anisotropic acceleration. 

In the present paper, we will firstly describe the formalism needed for discussing a 
general bigravity system. We will be interested in discussing cosmological solutions for 
such a system. As illustrative models we will use two types of potentials coupling the 
two metrics. One which resembles the Pauli-Fierz mass term in the linearized theory and 
one which is motivated by a higher dimensional brane construction. For these potentials 
we will be examining the most simple case, by imposing special symmetries, and without 
including matter. We introduce a description of the coupled cosmological evolution of the 
two metrics in terms of a "mechanical" model: two "relativistic particles" connected by a 
nonlinear "spring" . We will see that there is a generic period of acceleration of one or both of 
the metrics. For the Pauli-Fierz potential we discover two classes of accelerating solutions. 
One where anisotropies play a crucial role, and one where the cosmology is isotropic. On 
the other hand, for the brane potential our symmetry requirements only allow for solutions 
where anisotropies play an important role. In the final state the relative lapse between 
the two metrics tends to run away to infinity, when using the above illustrative coupling 
potentials. This run-away signals the breakdown of our effective theory. We discuss more 
general classes of potentials which naturally lead to a confinement of the relative lapse within 
a limited range. Such models lead to an interesting "locking" mechanism of the evolution of 
the two metrics. At the end, we briefly discuss the inclusion of matter in the above systems 
and propose the multigravity scenario as a candidate for a purely gravitational type of dark 
energy. 
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In the following we adopt the mostly plus metric signature (— , +, +, +) and use the fol- 
lowing definition for the Riemann tensor R K AMN = d M Ff N - d N T% M + T% n Tm H - T AM T% H . 
We use capital letters to label four dimensional spacetime coordinates and lower case let- 
ters to label three dimensional space coordinates. For the coordinate basis we use greek 
letters M, N, A, ... = 0,1,2,3 and [i,u,\,... = 1,2,3, and for the vierbien latin ones 
A, B,C, . . . = 0, 1, 2, 3 and a,b,c, . . . = 1, 2, 3. The reason for using, somewhat unconven- 
tionally, fj,,u,X,... for spatial indices will be explained below. 



2 General bigravity action 

A bigravity system is by definition the gravitational system of two coupled metrics which in 
linearized approximation is reduced to two spin-2 fields, one of which is massless while the 
other has a small massQ. The non-linear description of such a system can be achieved by 
considering the sum of the Einstein actions for two independent metrics (on the same man- 
ifold) and in addition an ultralocal potential term which couples the two metrics (see |[23|). 
The action is invariant only under the group of common spacetime diffeomorphisms. This 
guarantees the presence of one massless spin-2 excitation in the gravitational spectrum^]. 
Adding the coupling of matter fields {$1} to the first metric and of fields {$2} to the second 
one, we have the action: 

S = S G + S M (1) 



with 



Sg = ^J ^^ [gl] + 2^ / ^^fe] - V* J (9i92) 1/4 V (2) 



Sm = J v / =^TA(gi,{<fi}) + J v / =^A(g 2 , {$2}) (3) 

where V is a dimensionless scalar function of the relative metric gj~ 1 g2 and [ki] = [k 2 ] = 
M~ 2 while = M. Hence, the generic bigravity model has three dimensionful parameters, 
Ki, k 2 and fi. In the following we will factor out from the potential term the quantity (2/c) -1 , 
where R = + k 2 ) denotes the average gravitational coupling, and absorb the resulting 

4 Actually the second mass can also be zero and the coupling between the two metrics will emerge at 
the non-linear level. For more details see p4| . 

5 Depending on the parameters of the bigravity Lagrangian the massless graviton can be made arbitrarily 
weakly coupled, but it always will be in the spectrum for all finite values of these parameters B4j. 
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parameter of mass dimension 2 in V, i. e. we work with: 

Sg = ^~ j V=9iK[ Sl ] + -L / v^^fe] - i / ( 9l92 ) 1 / 4 V (4) 
2«i 7 2k 2 J 2k J 

where now [V] = M 2 . We have also the obvious freedom to add cosmological terms for 
the two metrics, but these can be absorbed in the potential V. In the following, we will 
consider that the potential does not have these cosmological terms, in order to isolate the 
physics of the coupling part of the above potential. 

It is convenient to introduce special frames uj a = ej^dx M with respect to which both 
metrics are diagonal, i.e.: 

3 



d*? = 2>2V) a (5) 

A=0 

<^ = £A2V) a (6) 



A=0 



where Aq an d Aq 2 '* are negative. We are then interested in the relative eigenvalues of the 
two metrics: 

x(2) 

The potential V, in this notation, is a function of the A^'s, or equivalently of the //^'s 
(see [El] for more details). The induced energy- momentum tensor from the coupling term 
of the two metrics is for each of the two metrics: 

rpMN _ ^ 5S m fQ ^ 



with S m = — i J ((?i5 , 2) 1//4 V". In the non-coordinate basis co> A , it has the simple form: 



= + (-sum) 



(10) 
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where <j\ = J2a Pa- This gives rise to effective energy densities p\ = — T^\ , p 2 = — T^q 

and effective pressures P° = To®, P 2 a = T^ 2 ' a a generated for the two metrics by the 
coupling potential (and both measured in their respective physical units). 

The only restriction we shall impose on the form of V is that at quadratic level it 
should reproduce the Pauli-Fierz mass term if the two metrics are expanded around a flat 
background, since this is the only ghost-free Lagrangian for a spin-2 field. Apart from 
this, the form of V is rather unconstrained. For example, a possible potential term was 
considered long ago in the "strong gravity" theory p3l and reads: 



V = m 2 (-^ 1 )«-V4 ( _^ ) -i/4 { tr [((g-i - g^g,) 2 ] - (trKfo 1 - gr^gi]) 2 } (11) 

with u+v = ~. We should point out here that this particular choice of potential is, contrary 
to the potentials we shall consider, asymmetric in the exchange of gi and g2- 

In the general case in four dimensions, the potential depends on o"i, a 2 , 03, 04, i.e. 
V = V(<7i, 02, 03, 04), where a n = J2aPa- l n the following we will focus on the class 
of potentials which depend only on <j\ and 2 , i.e.. V = V(ai,a 2 ). Note that for these 
potentials, (|) yields the following energy density and pressures for the first metric: 

Pl = ±e CT i/ 4 (\V - d ai V - 2/i d a2 V) (12) 

pa = _l e a 1 /4 Qy _ Q^y _ ^ Q^yj (13) 

while ( |I"0|) yields for the second metric: 

p 2 = \e~^ (\V + d ai V + 2/io d a2 V) (14) 
pa = _I e -«i/4 (|y + d ai V + 2fi a d a2 V) (15) 

This class of potentials is the minimal class which can reproduce the Pauli-Fierz mass 
term in the limit where g\j N and g 2 MN are expanded around the same flat metric r] MN . The 
particular combination of 01 and a 2 which accomplishes this is: 



rn 
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V{a 1 ,o 2 ) = ^{o 2 -ol) (16) 
Indeed, expanding the two metrics around a common flat background as g\ [N = t/mn + 



y/2i^h\ IN and g\ IN = r] MN + ^2K 2 h\ 1N , and introducing x x = . and x 2 



K1+K2 A/ K1+K2 



K2 



satisfying x\ + x\ = 1, we find that the combination h° MN = X2h l MN + x\h 2 MN is massless 
while hZ N = x ih\i N - x 2 h 2 MN = , 1 \ g\ 1N ~ g 2 MN ) has a Pauli-Fierz mass term: 

- Yn {9l92)1/4V = " kl) (17) 

The above potential is the simplest choice that one could think about and hence we will 
mainly focus on it in the following. We shall call the above potential the "Pauli-Fierz 
potential" and examine the physics when the quantities <j\ and a 2 are not necessarily 
small. 

We shall also consider the potential that arises as the effective four dimensional descrip- 
tion of the brane motivated bigravity scenario JTJ. This potential has the form [ 23J : 



V = m 2 [cosh B — cosh A] , (It 



with A = — and B = —^=\ a 2 — — (19) 
4 2V3\ 4 

The potential V in this case is a function solely of ui and a 2 - In addition, it is symmetric 
under the exchange of gi and g2 (or equivalently it is an even function of the /m's). For 
small o\ and a 2 we get exactly the Pauli-Fierz potential (|16D , with the change of notation: 



m PF (20) 



The discussion in |M] of the Arnowitt-Deser-Misner Hamiltonian approach to bigravity 
has shown that a crucial feature of the potential ^ A (g]~ 1 g2) is its ability (or lack of ability) 
to "confine", within a limited range, the variation of the relative lapse n = (N 2 /N 1 ) 1 ^ 2 , by 
means of its (algebraic) equation of motion. We have found that the confining capabilities 
of both the simple "Pauli-Fierz" potential (|T6|) and the "brane-motivated" one ([18]) are 
too weak to prevent the relative lapse (appearing as 7 = 21ogn in our notation below) 
from running away to infinity in a finite (proper) time in generic solutions. This led us to 
introduce and study modified versions of the above two potentials. For instance, instead of 
the purely "quadratic" (in the /i^'s) Pauli-Fierz potential (|i~6|) , we shall consider a Pauli- 
Fierz potential augmented by "quartic" terms, say (with the notation m = \/3 tupf)'- 

2 

V(a 1 ,a 2 ) = —(a 2 -a 2 1 +Xa 2 2 ) (21) 
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where A is a positive parameter of order unity. 

In the above cases, even though the potentials are symmetric in g x and g2, they can 
couple asymmetrically to matter if k\ and K2 differ. For the sake of simplicity in the 
following we constrain our study in the case where K\ = k 2 , where we obtain what is known 
as symmetric bigravity. Additionally, we normalize K\ = K2 = 1/2 and leave the only 
dimensionful parameter of the effective action to be m = a/3 mpp. 

3 Cosmological ansatze 

For the above choices of potential we wish to study the cosmological evolution of the 
bigravity system. A first assumption that we make is that the two metrics depend only on 
time. Then the line elements will have the form: 

ds\= -e 2 ^{dx ) 2 + xlvidx^ + b^dx°){dx u + b u lC lx ) (22) 
ds\= -e 2 ^{dx ) 2 + xlu(dx^ + b^dx°)(dx v + b u 2 dx°) (23) 

where all the functions indicated are only time-dependent. The residual symmetry of such 
cosmological metrics are the following common diffeomorphisms: 

x°= f(x ') (24) 
x »= x^' + ^(x ') (25) 

where f(x ') and ^(x ') are arbitrary functions of time. Under such common diffeomor- 
phisms the various fields transform as: 

P 2 7 i _, p 2 7l f/2 2 72 __. p 2 7 2 f/2 1 1 2 2 

° ° J 1 ° ° J 1 Afiv Afiv 1 Afiv Afiv 1 

bi^Kf + (ey , b^btr + iey (26) 

so that the quantities 7 = 72 — 7i, xL, xjL an d ^ M = e_ ^(^2 — ^i)> where 7 = 71 + 72 , are 
invariant. If we where to consider only the Einstein terms in the action (|2|) together with 
the matter terms (|3]), b^ and b^ would only enter (modulus surface terms) as Lagrange 
multipliers in front of the momentum constraints. As these constraints identically vanish 
within the simple ( "Bianchi Type I" ) cosmological solutions that we consider, we conclude 
that 6^ and b% only enter the action through the potential V. Moreover, because of the 
transformation properties exhibited in (j2"6|) , the invariant V can only be a function of 7 = 



7 



72~ 7i) xL/5 Xuw an d the combination 6 M = e~^ {b^ — b^) with 7 = 71 ^ 72 as above. Taking into 
account the fact that detg\ IN = — e 271 detxL, (similarly for the second metric), we see that 

the shift vectors enter the action only through e^(detx^) 1 ^ 4 (detx^) 1//4 ^(7, b^, xL/> X/D- 
As V is a scalar, it can involve the vector 6 M only through some scalar combinations such 
as h^b v , k^ UKX b^b u b K b x , . . ., where h^, k^ UKX , ... are made from x 1 ^, X%, ancl 7- As V 
is by assumption a smooth function of gj~ 1 g2, it will be a smooth function of the scalars 
made with b^, so that we can write: 

V(b») =V + h^b v + k^ VKX Wb K b x + 0(b 6 ) (27) 

where V , h^ u , k^ UKX , depend only on x\ v i X% an d 7- The equation of motion for bi, b% is 
simply that 6 M extremize V^), i.e. J^- = 0. We see immediately from ( |2"T| ) that ¥ L = 
is always a solution of this equation. In this work, we shall only consider this universal 
"perturbative" solution. Note, however, that, similarly to what happens in Landau's theory 
of magnetic phase transitions, there might, for some potentials, exist also some "nonper- 
turbative" (or "symmetry breaking") solutions with U L 7^ 0. In Appendix B we study the 
equation of motion of 6 M by slightly different approach, and find that "nonperturbative" 
solutions can only exist if the potential V(0"i,o"2) is such that d a2 V can vanish (which is 
the case neither for ([16]) nor for (|18D). 

Once we know that b^ = 0, i.e. b 1 ^ = b%, we can further use the ^(x ') gauge freedom 
to set both bi and b% to zero. Then the metrics become: 

ds\= -e 2 ^{dx ) 2 +xU^dx u (28) 
ds 2 2 = -e 2 ^{dx ) 2 +xl„dx»dx u (29) 

Note that we still have the gauge freedom f(x 0r ) which changes 71 and 72 but leaves 72 — 71 
invariant. 



4 Cosmological dynamics of a bigravity system 

The action describing the dynamics of xL,, x\ v (obtained by inserting the ansatze (|28|), 
fl29|) into the action (f|); without the matter terms for simplicity) reads S — J dx°C with: 

Co = e-^/dfri, *i) + e- 72 /C 2 (x 2 , X2) ~ e^+^V 12 (7 2 ~ 7i, Xi^) (30) 
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where the kinetic terms are (using matrix notation for x\v> xfm)'- 

MXi, Xi) = ^(detxi,) 1/2 Hxl'Xx) 2 ~ (Hxl'Xi)) 2 } (31) 
with a corresponding definition for /C 2 (% 2 , X2) where xL — * xjL, an d: 

V 12 ( 72 - 7i, Xi^) = (detxi,) 1/4 (det X ^) 1/ V( 72 - 71, X^Xi) (32) 

If we replace the two variables 71, 72 by the equivalent combinations 7 = 5(71+72), 
7 = 72 - 7i (so that 71 = 7 — ^7, 72 = 7 + it is easily seen that 7 is a gauge variable 
whose equation of motion gives the usual zero-energy (Hamiltonian) constraint: 

B = e^fcxfo, xi) + e- 7/2 /C 2 ( %2 , X 2 ) + V 12 ( 7 , X^Xa) = (33) 

After imposition of this (conserved) zero-energy constraint, we can use the f{x 0> ) gauge 
freedom (pf ) to set 7 to zero say By contrast, 7 is a dynamical variable whose equation of 
motion is algebraic, because 7 has no kinetic terms of its own. The equation of motion of 
7 is that £0(7) should be extremized, i.e. ^ a = 0: 

j- [e^ixt, Xi) + e" 7/2 /C 2 (x 2 , X 2 ) ~ Va(7, X^Xa)] = (34) 

We shall assume (see below) that the evolution remains in a domain which confines 7 to 
a limited range, i.e. that the shape of the potential V(7) and the values of the kinetic 
terms /Ci, /C 2 , are such that there exists a bounded solution 7 of equation (|34D which can 
be continuously followed during the time evolution. Then the main problem is to discuss 
the coupled dynamics of the matrices and xjL, obtained (say in the gauge 7 = 0) from 
the Lagrangian: 

C = e^ 2 /d(Xi, Xi) + e~ 7/2 /C 2 (x 2 , X 2 ) - Vi 2 ( 7 , XT 1 *!) (35) 

For obtaining the equations of motion of x\v and x\v> one can equivalently either treat 7 as 
an independent variable, which will be later replaced by the solution of equation ([51]), or as a 
function of Xi-> Xi, X 2 > X 2 defined by equation (|34"|). The latter way of viewing 7 transforms 
(|35f) in a very complicated and very nonlinear Lagrangian £ re duced(Xi, Xi> X2; X 2 )- For 
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simplicity, using our assumption that 7 remains bounded, i.e. 0(1), we shall view 7 (until 
it is replaced in the final equations of motion) as a "given" function of time. 

In this simplified view, the dynamics of x]iui. x °) an d xfiu( x °) can be interpreted as a 
mechanical model of two "particles" (with 7-dependent, i.e. time-dependent masses), liv- 
ing in a six- dimensional pseudo-Riemannian space, and connected by some type of (time- 
dependent) nonlinear "spring". Indeed, each symmetric 3x3 matrix i — 1,2 has 
six independent components, and the kinetic terms ([H]) define a certain (curved) Rie- 
mannian metric, with signature (—,+, +,+,+, +). The (7-dependent) nonlinear "spring" 
is defined by the potential V^ix^X-i) i n equation fl32|) . It is interesting to note that, 
in the case of the brane potential, the potential Vi^Xi" 1 ^) happens to be precisely 
the squared geodesic distance between Xi an d X2-> as defined by the Riemannian metric 
dx 2 = A^(x, dx) = K-(XiX)dt 2 P^fl . Therefore, in this case, modulo the 7-dependent 
modulations, the cosmological evolution system can be elegantly viewed as the problem 
of two "particles" in a pseudo-Riemannian space, connected by a harmonic "spring" as 
Vi 2 (xr 1 X2) °c (distance) 2 . 

In spite of the possibility of such an elegant formulation, the actual dynamics of the 
coupled metrics x\v an d x\ v * s extremely complicated. Therefore, as a first cut towards 
understanding the main qualitative features of this dynamics, we shall henceforth further 
specialize the class of solutions we consider by focusing on the metrics which can be simul- 
taneously diagonalized. It is then convenient to parametrize the diagonal components of 
xliv as x\v = ^ a S ( e2aM )' an d those of x\v as X% = diag (e 2/3M ). In other words, the two 
metrics read: 

3 

ds\= -e 2 ^{dx ) 2 + ^e 2a "(^) 2 (36) 

3 

ds\= -e 2 ^{dx ) 2 +^e 2/3 "(^) 2 (37) 



Note that this class of metrics is more general than the "bi-Friedmann" metrics where 
and f3^ would both be taken to be (|, |, |) and (|, |, |) where a = ^] ct M and (3 = ^ (3^. 

The consideration of unequal a M , (3^, [1 = 1,2, 3, means that we are considering anisotropic 
metrics. As we shall see, anisotropies can play a crucial role in the dynamics of "bi- 
cosmology" . 
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For the restricted metrics (|36|) , (|37D the Lagrangian fl3~Cf ) simplifies to: 



C = e a ~^ G^a u + e^ 2 - e {a+ ^ + ^ )/2 V( l2 - 7l , /3" - a' 4 ) (38) 

where the contractions are made with the internal metric: 



sx«" = ^) J - fey) 



(39) 



Note that the inverse metric reads: 



(40) 



where the factor 1/2 would be 1/(D — 2) in spacetime dimension D. We note here that 
the above metric has a Lorentzian signature (—,+,+) and thus we can think of a M and 
/5 M as the worldlines of "particles" in a three-dimensional spacetime (instead of the above 
six-dimensional spacetime). The use of greek indices \x = 1, 2, 3 for labeling the three inde- 
pendent scale factors a M = e" M was chosen to emphasize the Lorenzian nature of the metric 
Q^yda^da 11 underlying the kinetic terms in (|38"D . Actually, the "Poincare" symmetry in 
field space (under transformations preserving the scalar product defined by Q^v) of the La- 
grangian (^) is broken by two sets of terms. Firstly, the separate appearance of a = ^ a M 
and (3 = 5^/?^ in the kinetic terms introduces a preferred covector = (1, 1, 1) such that 
a = n^a^ and (3 = n^^ 1 Q This covector is timelike because n 2 = Q^n^n^ = —3/2. 
Actually, we see that, in conformity with the metric dx 2 = K>{x, dx) introduced above in 
the space of metrics x M „, the Lagrangian ([3~8~D features the curved metric e a Q^da^da u . We 
found, however, convenient to think in terms of the (conformally related) flat metric Q^ v . 
Secondly, the potential V(a 1 ,a 2 ) depends on: 

o x = 2( 72 - 7l + p - a) (41) 
a 2 = 4 ( 72 - 7l ) 2 + 4E(^-^) 2 (42) 



6 Note that the notations a = n^a^, j3 = and later a = n^a^, S — n^S^ correspond to some linear 

projection of the "spacetime" vectors on a certain timelike direction. They could have been denoted ao, 
Po, etc.. They should not be confused with the magnitude of the corresponding vectors Q^a^a 11 which we 
shall never encounter below. 
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which also breaks the field space "Poincare" symmetry. It is crucial to note that, as in- 
dicated in (^), V depends only on the differences 72 — 71 and [3^ — ct M . This suggests to 
introduce the sums and differences of the above particle coordinates to be: 



(43) 
(44) 



Using the f(x ') diffeomorphisms fl24]) , we can as said above fix the gauge ambiguity by 
imposing 7 = §(71 +72) = 0, i.e. by setting 72 = —71 = 7/2. In this gauge we are using as 
time parameter the "average proper time" t, with dt = e^dx . Then after further defining 
a = ^ cr M , 5 = ^ 8^, we obtain the action S — J dtC t with: 



C, 



D a/2 



,a/2 



I cosh f 1 ^) + - sinh ) - 



(45) 



where the contractions are made with the internal metric (0). When using such a gauge 
fixed Lagrangian one must remember to impose separately the constraint coming from the 
variation of the average lapse function 7 = |(7i + 72). As said above, this constraint is 
the zero-energy condition, i.e. £ = K + V = where K is the total kinetic energy. As 
said above, among the equations of motion, two play a special role in that they do not 
involve second order derivatives. These are the equations of motion obtained by varying 
the two lapse functions 71 and 72 (which are essentially the Friedmann equations for the 
two metrics). They read (in the 7 = gauge): 



V _ 9V ) -( 7 -5)/2 
2 d-y I 



(46) 
(47) 



which, using the definitions <^), ([UJ), is equivalent to the more suggestive expressions: 

da^ da 



e^a^a^ + p\ = 



dt\ dt\ 

dp d(3^ 
dti dti 



+ Pi = 



+ p 2 = 



(48) 
(49) 
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where dt\ = e J ' 2 dt = e^dx and dti = e'^dt = e l2 dx° are the separate proper time 
coordinates for the two metrics. 

Alternatively, working in terms of the combinations 7 = -(71 +72) and 7 = 72 — 71, the 
equations of motion obtained by varying respectively 7 and 7 read: 

e- a/2 £ = ~ cosh [p^J + 6%) - sinh f 2 ^) + V(j, 5 M ) = (50) 



V = 2 smh (V j ^ + ^ - cosh (V j °* 6 » — 



The above two sets of equations, (0), ( ^9|) versus (|50"D, (|51|), can be easily checked to 
be totally equivalent. Note that while equation floTf ) is really a first class constraint, which 
is preserved by the evolution equations of motion, and only needs to be enforced at some 
initial time, the equation (]51f) is an algebraic equation which determines 7 as a function of 

a^, d M , /3 M , /3 M or equivalently a M , d M , <5 M . As said above, unless explicitly mentioned 
otherwise (see below), we assume that we are in a regime where 7 can be continuously 

solved in terms of cr^, (T M , <5 M , 5^. Then the main evolution system is obtained by varying 
a^, or equivalently cr^, in the Lagrangian (f45|). See Appendix A for the explicit 
form of these evolution equations. 

As we said above, the Lagrangian involves, both in the kinetic terms and in the potential, 
a special covector n M = (1,1,1). This makes it useful to decompose the above functions 
into timelike and spacelike components as: 

5" = -§n"<5 + (52) 
a* = ~n"a + (53) 

with the property that Wn^ = S 1 + S 2 + S 3 = and A^ = A 1 + A 2 + A 3 = 
0. In the definitions (|52]) , ( |53"1) there appear the contravariant components of n^, i.e. 



n fj. = g^ nv = —-(1,1,1) and we recall that n^n^ = — |. We can easily check that 



the longitudinal components a and 5 appearing in (|52"D, ( p3| ) are simply cr = cr^n^ and 
5 = 5^7^. It is helpful to keep in mind the geometrical configuration of Fig.|l] which 
displays the various bimetrical parameters which are relevant for defining the dynamics of 
bigravity. The elements of this geometrical configuration which play a dominant role in 
our discussion are: (i) the timelike versus spacelike character of the "velocity vectors" d' 4 , 
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or 



Figure 1: General motion of the two "particle" worldlines and (3^ in the internal 
Lorentzian field space. The middle worldline is the "center of mass" of the system 
The decomposition of the "center of mass" velocity a fl /2 and of the worldline separation 
5^ in timelike and spacelike components is shown. These projections are performed with 
respect to the special timelike vector n M . 



/3 M (with respect to the light cone defined by Q^), and (ii) the timelike versus spacelike 
character of the "separation vector" 5 M = (3^ — a M . 

We can define now the fields: 

3 

= 2< 
V3. 



A i 



A_ 



-(A 2 + A 1 ) 
-(A 2 - a 1 ; 



(54) 
(55) 



and respectively for E + , S_. Then the Lagrangian can be expressed as a function of the 
timelike fields a, 5 and the pairs of spacelike ones A±, E± as: 



,a/2 



- cosh f - — - 
3 V 2 



-a 2 - 5 2 + Ai + A 2 + Ei + E 2 



~ sinh f^A ) + A+ E + + A_E_) - V(j, 5, A +! A_) 



(56) 



As indicated, a generic potential can only be a function of 7 and the 5 M 's, i.e. of 7, 5, 
A + and A_, hence we see that E + and E_, which correspond to the "spatial" location of 
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the center of mass of the two particles a^, (3^, are ignorable coordinates. In fact, in the case 
of our class of potentials V(a±, <7 2 ), the dependence of V is even more restricted because: 



a 1= 2(7 + 5) 



(57) 
(58) 



where r 2 = |A^A^ = A 2 + A 2 , so that the potential depends only on the magnitude of 
the "spatial distance" separating the two particles. It is convenient to go to the Routhian 
formalism and replace the S + and XL by their conserved canonical momenta. The Routh 
functional is then: 

-TZ t = - £-Pe_ + C t 

3 (p| + +p|_ 



= W 2 



+ - cosh f - — 
3 V 2 



A 2 + A 2 



TV"- 2 -* 



where the conserved momenta 



cosh 2 (^) 
are: 



dCt 2 CT/2 



cosh ( — — -J S-t — sinh 



A. 



(60) 



For simplicity we shall focus on the case where p-£ + = = 0, i.e. when the Routhian 
has no terms linear to A + , A_. This is the special case where the "total momentum" of the 
two-particle system is directed along the "vertical" time axis n M . As the formal "Poincare" 
invariance of the dynamics is broken, note that going to such a "center of mass" frame 
is a real restriction on the class of solutions that we examine. Then there exists another 
conserved quantity because of the rotational symmetry of the system with respect to n^. 

This can be seen by defining: 



A i = r sin 9 



A_ = r cos 9 



(61) 



Then the Routhian is written as: 



7?, r'U-m^J' 



-.2 _|_ r 2g2 

cosh 2 (^) 



a 



- sinh 

3 



7 



&5-V{-f,5, r 2 ) 



(62) 
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The potential depends only in r, thus there exists a conserved "angular momentum" pg 
given by: 

&R T 2 r 2 6 



Pe 



86 "3 cosh (2=2) 



,<r/2 



(63) 



In analogy with our previous treatment, we can define a new Routhian eliminating the new 
ignorable coordinate: 



Tit 




+ ~sinh (l-l)a5-V( 1 ,5,r 2 ) 



^2 rz 

a — o 



cosh 2 pz£) 



(64) 



In the most simple case we can have also pg = 0, where the motion of the two particles 
is planar. In this case the variable r (which represents say, A_, if 9 = const. = 0) can be 
considered as varying on the full real line, from — oo to +oo. We will mainly discuss this 
particular case and from now on, we will indicate by £ the simplified Routhian: 



£ ,-w4ico sh (X^) 



,j. 2 



cosh 2 ( V) 



& 2 -S 2 



(65) 



Let us recall that the Lagrangian (|65D corresponds to the restricted configuration where 
both particle worldlines lie in the same timelike plane. The logarithms of the relative 
eigenvalues of the two metrics in this special limit are: 



IM) = 27 



\i\ = -(6 - VSr) 



(66) 



Even though up to now we have restricted the conserved quantities that appear in the 
system to have particular values (zero), the dynamics of the above system is still very rich 
and difficult to analyse. To get an insight in the behaviour of the solutions of the system, 
we will study two extreme cases. In the one case the two world lines and (3^ are mostly 
spacelike separated, and in the second mostly timelike separated. An extreme version of 
the first case is realized by setting 5 = and having r ^> 1. Similarly, an extreme version of 
the second case is realized in the opposite case where r = and 8 ^> 1. Moreover, we shall 
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further restrict the scope of our analysis by requiring that the velocity vectors be, at least 
initially, both timelike (and future directed). We will first study the case of the Pauli-Fierz 
potential and then indicate the changes when considering the brane motivated one. 



5 The Pauli-Fierz potential 

The Pauli-Fierz potential ([16]) with the previous notation has the simple form: 

2 

HI 

V( 1 ,5,r) = — (r 2 - 5 2 - 35 1 ) (67) 

Let us mention the values of the energy density, the pressures as well as the equations 
of state for the different spatial directions in this case for the first metric: 

2 

Pl = ^-e 2 ^ (r 2 - 5 2 - 3 7 <5 + 66) (68) 
18 

2 

P«= -^-e 2 ^ (r 2 -5 2 -3 7 5 + 67 + 45-2v / 3{-l,l,0}r) (69) 

The corresponding quantities for the second metric are respectively: 

p 2 = ^e-^(r 2 - 5 2 - 3 7 5 - 65) (71) 

2 

P« = -^e" 2 ^ (V 2 - S 2 - 3 7 5 - 67 - 45 + 2>/3{-l, 1, 0}r) (72) 
18 V / 

a_ 1 _ /-3o^ 2 V3{-u M o }f 

P2 ^ — — 3 7 o — 60 

Let us now specialize in the two aforementioned limits. 

5.1 Spacelike separated worldlines 

We consider first the purest spacelike separation corresponding to the limit 5 = 0, where 
the two worldlines are symmetric with respect to the motion of their "center of mass" (see 
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Figure 2: Motion of the two "particle" worldlines a M and (5^ for the exact spacelike worldline 
separation limit. The middle worldline is the "center of mass" of the system and is 
directed along the special vector vV\ The motion of the two "particles" is symmetric with 
respect to the motion of their "center of mass" . 



Fig.0). It is easily seen that S = is always a solution, if and only if one has also 7 = 0. 
The 7 and 5 equations of motion are then satisfied independently of the evolution of the 
other degrees of freedom. The potential is simply: 

2 

77? 

V{r) = — r 2 (74) 
and the Lagrangian (^) simplifies to the rather friendly form: 

C = e°' 2 (-r 2 - -a 2 - —r 2 ] (75) 
V 3 3 9 J 

The above Lagrangian nicely illustrates the view of the dynamics of the system as that of 
two particles connected by a spring (with an energy proportional to r 2 , i.e. to the square 
of the spatial distance separating the two particles). 

The system of equations of motion for the fields a and r is: 



2 /. a 



3 



a 



V(r) (76) 
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with the Hamiltonian constraint: 

-(r 2 - a 2 ) + V(r) = 
3 



(78) 



It is easy to check that by differentiating the Hamiltonian constraint and using one of the 
equations of motion we can obtain the remaining equation of motion. The structure of the 
above system is similar to the one of a scalar field in an isotropic spatially flat universe. 
The scalar field in the above case is r while a is the logarithm of the scale factor. Thus, if 
the slow roll conditions \d r V jV^ <C 1, |Qj?V/V| <C 1 hold in the above system, there will be 
at least one inflationary domain. In the particular potential at hand these two conditions 
can be satisfied when r>l. 

In this slow- roll region we can find solutions of the fl76f), (|77|) , (f78|) by omitting the 
second order derivatives in the equations of motion and the f 2 term in the Hamiltonian 
constraint. The equations can then be easily integrated to: 



r 2m + 



a = a + 



r t 



m j.2 
V3 l 



(79) 
(80) 



where the relation & = ^=r holds. The logarithms of the scale factors in this case for 

the three spatial directions are (after a rescaling of the spatial coordinates to absorb an 
arbitrary integration constant): 



ct I r 
3 ^ V3 

- { £ _ J_ V and 8^ 
2 1 3 y/3 ' 



1 

2< 



3 V3 
ct _i_ r 
3 v/3 



Thus, from the above formulas we see that the two metrics will be exponentially inflating 
until the time when r ~ 0(1). Note that the solution is written in proper time for both 
metrics since 7 = 0. The anisotropic Hubble parameters for the two metrics are defined as 

= d M and = (3^ and are explicitly, along the three spatial directions: 



Hi 



6t 



m 



r 
2^3 



> and H% = 3» 



m 
3" 



-j— 4. 1 

2V3 + 

-L 1 

2V3 

r 
2V3 



(82) 
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Then, the mean Hubble parameter will be the same for the two metrics and equal to: 



< H X >=< H 2 >= ± J>f = \H H Z 



a 
6 



(83) 



Thus, we can quantify the anisotropy of each metric by the introduction of the parameters: 



\ 



3 < Hi > 2 



A 9 = 



^ 3<H 2 > 2 



i4) 



which in the particular case at hand are: 



A 1 = A 2 = V2 



r 
a 



2V2 



15) 



Hence, the anisotropy in each metric increases during inflation and becomes 0(1) when 
the slow-roll condition is violated. This is an interesting prediction of this type of bigravity 
"inflation" which drastically differs from the scalar-inflaton-driven inflation which tends to 
wash out any initial anisotropies. 

The equations of state for the first metric and for the three different spatial directions 
during the inflationary era, and their limit for r ^> 1 can be read from (|70|): 

pa ^ 

wl = — = -1 + 2v / 3(-l,l,0)- -> -1 (86) 
pi r 

and for the second metric w 2 are the same with a flip in the sign in the second addendum. 
Note the simple relation that the equations of state for the two metrics satisfy: 

wl + w « = -2 (87) 



When r ~ 0(1) our approximation breaks down and r will start to oscillate around 
zero. In this limit a grows logarithmically. The anisotropy parameter in both metrics is 
maximal (and 0(1)) near the zeros of r, and zero at the extrema of r. The evolution of the 
system during the inflationary and the oscillatory region is illustrated in Fig.|3|. However, as 
we are going to see, the later oscillatory region will cease to exist in the presence of a slight 
initial perturbation of 5 and the relative lapse 7 will undergo (for the chosen potential) a 
run away towards infinity before the first zero of r is reached. 
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Figure 3: Numerical simulation of the motion of the two "particle" worldlines a 11 and (3^ 
for the exact spacelike worldline separation limit and for the Pauli-Fierz potential. The 
horizontal axis is the distance r/2 from the "center of mass", while the vertical one is a and 
/3 respectively for the two "particles" . Initially the worldvolume of both metrics is inflating 
until the worldline separation becomes 0(1). Then the two "particles" start oscillating 
with respect of their "center of mass". The oscillatory regime cannot be reached for the 
Pauli-Fierz potential if the timelike separation component 5 is excited. 



Let us indeed consider the effect of perturbations of 5 away from 8 
7 = 5 + 2e, then the perturbation of the Routhian in quadratic order is: 



0. If we let 



5C = e 



a/2 



2 -ea'b- lm 2 (25 2 + 35e) 
3 9 



The constraint ^ = can be used to solve for e. If we substitute it back to the Routhian 
we find: 



3 \ (J -\- T 



5 2 - —e^ 2 
3 



'4 2m 2 -a 2 fa 
- H 2 

3 a 2 + r 2 \& 2 + r 2 



(89) 



B 
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The quantities A and B during the inflationary period and their limits for r ^> 1 are: 



B = ^ 2 



r 2 - 10 







3 r 2 + 4 


3 








) 8r 2 






(r 2 + 4) 2 _ 
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(90) 
(91) 



Thus, the perturbation of the action simplifies (during slow-roll) to: 

5S = J dte^ 2 (-i<5 2 



(92) 



From the extremization of this action we get the following motion for 5: 



S = e^ 2 



C x + C 2 Erf 



2V2 



(93) 



where Erf(x) = ^= J* e y dy is the error function. Thus we see that although the field 

S has an unstable looking action (after taking into account the ghost-like kinetic term of 
S its correct-looking mass term corresponds in fact to a tachyonic instability), the time 
dependence of o tames and actually damps the evolution of 5. This is confirmed also by 
numerically integrating the full equations of motion for a small 5 perturbation. 

However, in the region where r ~ 0(1), where the result of the above perturbation 
analysis breaks down, we see numerically that 7 starts growing fast and finally, always 
before r reaches its first zero, it runs away to +00. In this limit we see that a^a^ — > 0~ 

while f3^(3^ — > +00. Note that remains always timelike (a^a^ < 0), but (3^ starts timelike 

(frfin < 0), goes through its "light cone" (/^/^ = 0), and ends up spacelike {$ fl ^ > 0). 
To understand what happens in the limit 7 — > +00 we have examined numerically three 
invariant quantities for the two metrics. The first is the proper energy densities: 



pi = —e'a^a^ 



P2 = 



1-1 w 



(94) 



which are proportional to the Gq component of the Einstein tensor. The limiting behaviour 

of a^a^ and (3^/3^ conspire in such a way with the runaway of 7, that the above energy 
densities tend to finite values, positive for p 1 and negative for p 2 . The second invariant 
quantity is the anisotropies of the two metrics. For the first one we observe that after 
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having reached a maximum value, it starts decreasing, and for the second one we observe 
that close to the point where 7 diverges, it has a decreasing behaviour as well. The third 
set of invariant quantities that one can study is the equations of state. For the first metric 
all three equations of state have the limit — > — 00 and for the second metric they tend 
to negative but finite values. From all the above we conclude that the final state of the 
system as measured by each metric is not singular even though the relative lapse 7 runs 
away to +00. 

This running away, however, signals that our effective theory description breaks down 
and that the limit we are discussing is governed by different dynamics. Indeed, the mass 
scale corresponding to m 2 becomes m\ ~ m 2 e 7//2 when viewed in the first metric and 
mj ~ m 2 e _7//2 when viewed from the second one (see e.g. equations (|68|), (|7l|)). Therefore, 



when 7 gets too large, this initially "light" scale becomes "heavy" on the first "brane"Qand 
"ultralight" on the second. This means, for instance, that the tower of heavier graviton 
masses above m\ (that are treated as infinitely heavier and have been truncated away in 
deriving our effective four-dimensional action) might become light and should be taken into 
account. The run away of 7 then signals the necessity to shift to a new basic Lagrangian. 

A simple reason of the run away of 7 in the above case can be seen from the structure 
of the original Lagrangian (|45|): 



2 

m 



C t = e a+ ^ 2 + e^ 2 - ^ 2 (r 2 - 5 2 - 3<5 7 ) (95) 

9 

Seen as a function of 7 the above Lagrangian has three important parts, namely an 
increasing exponential in front of the kinetic term of the one "particle" , a falling exponential 
in front of the kinetic term of the second "particle" and a linear contribution in the potential. 
If the kinetic terms of the two particles have the same sign (i.e. if they are both timelike 
or both spacelike), there will always be an extremum of the action for a finite 7. However, 
when they start taking opposite signs the extremum can easily run away to infinity and 
disappear (except when 5 has the good sign and is large enough to confine 7 in some 
bounded interval). 

As said in Section 2 this led us to consider modified potentials, with improved "confin- 
ing" properties for 7. For instance, an improved version of the naive Pauli-Fierz potential 



7 We do not have in mind any brane configuration that can lead to a pure Pauli-Fierz potential in its 
effective four dimensional action. However, we can still think of "branes" being weakly coupled worlds 



(see 1 24 ) in a higher dimensional setup having this particular four dimensional effective description. 
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fll6f) is the potential ( plf ) which includes a quartic contribution (note that one generically 
expects the presence of such quartic contributions in nonlinear potentials V(iia)): 

2 

V(a ll a 2 ) = —(a 2 -a 2 1 +Xa 2 2 ) (96) 

Now the potential fl6"7|) will include in addition to the linear in 7 term, a term which is 
quadratic and a term which is quartic in 7. These terms help to ensure that there exists 
an extremum for finite 7 and thus that the "particle" velocities can go through the light 
cone without any pathologies. Note, however, that the existence of a solution in 7 is not 
guaranteed for arbitrary evolutions. Indeed, the potential ^(7), which tends, by itself, to 
confine 7 near 7 = 0, is compounded by the effect of the kinetic terms in equation (|D~5|) 
which become actually dominant for large values of 7. When a^a^ and (3^(3^ have opposite 
signs these kinetic-related exponential potentials might destabilize the "local" confining 
ability of V ('-/). The only case where one would be guaranteed to always have a solution 
in 7 is the case where a confining V(7) would dominate over e' 7 '/ 2 as I7I — > 00. Anyway, 
we find that the slight deformation of the potential brought by the quartic addition Xa 2 , is 
sufficient (for the solutions we explored) to prevent the runaway of 7 and the breakdown of 
the effective theory in the pure Pauli-Fierz potential. Consequently, the oscillatory region 
will be present for ( ]9^| ) even if 5 is slightly perturbed. We have confirmed the above 
behaviour numerically and further observed that the maxima of the anisotropy parameter 
during successive oscillations are of decreasing amplitude. 



5.2 Timelike separated worldlines 

An extreme example of timelike worldline separation is the limit r = 0, where the two 
worldlines are collinear (see Fig.|]). One checks that if initially r = f = 0, it remains zero 
all over the evolution. The potential becomes: 



V( 7} 6) = -^{6 2 + 35 1 ) 



(97) 



Then the Lagrangian (pSI) simplifies to: 



,a/2 



- - cosh I - — - 
3 V 2 



o 2 + S A ) + - sinh 



7 



&S + ^- (6 2 + 35 7 ) 



(98) 
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Figure 4: Motion of the two "particle" worldlines a M and (3^ for the exact timelike worldline 
separation limit. The two worldlines are collinear and the evolution of the two metrics 
isotropic. 



Since setting the spacelike separation between the worldlines to zero actually makes 
the two metrics isotropic, we can express the above Lagrangian in the following equivalent 
form: 

C = -\e a+ ^ 2 a 2 - le^ 2 $ 2 + e^ 2 ^ (5 2 + 35 7 ) (99) 
So 9 

From the above it is worth noting that the S 2 term of the potential has tachyonic sign 
while both a and S (or a and j3) have ghostlike kinetic terms. This means that the potential 
V has again the "good" relative sign allowing stable motion in such a timelike configuration. 
In other words, we have now two "particles" which are vertically separated, and which are 
still connected by a confining "spring" (with potential quadratic in the vertical separation). 
We expect to have solutions where the particles will chase each other, overtake each other, 
etc. and go through some type of timelike oscillatory motion. This expected behaviour is, 
however, not realized in the case of the original Pauli-Fierz potential (|i~6|) because of the 
instability of the potential to confine the motion of 7. On the other hand, we have found 
that the expected behaviour is obtained when using the modified potential (|2T|). Let us 
start by briefly describing what happens in the case of the original Pauli-Fierz potential 

(HD- 

The Friedmann equations (f48D, ([49]) in this case demand the positivity of the energy 
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densities pi, p%. They read for this particular case: 




2 

m 



5(5 + 3(7 




(100) 
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2 

m 



5(5 + 3(7 



2))e 



2 



(101) 



18 



+ 



From these, we can easily see that there are two allowed regions in the (5, 7) plane, the 



allowed regions are away from the usual "perturbative" domain I7I (corresponding to 



The solutions fall into two categories. The one which gives accelerated expansion in one 
of the metrics and the ones that give decelerated contraction on the same metric. Since 
the solutions in the above-mentioned regions can be converted to one another through a 
time reversal or exchange of a and (3, it suffices to examine the first region only and in 
particular the expanding solutions. The motion of the system for expanding metrics in this 
case has the following behavior in the (5, 7) plane: initially 7 increases with time, while 5 is 
decreasing. If we evolve the solution back in time, we hit a cosmological initial singularity, 
where 7 ^ — 00, 5 ^ 00. If we evolve forward in time, we find that, after 5 has decreased 
to the region where 5 ~ 0(1), 7 reaches a maximum value and starts decreasing, while 5 
keeps decreasing until 5 — > 0. At this latter limit 7 — > — 00. In both regions a is increasing 
in such a way that (3 is very slowly varying being in very good approximation constant (i.e. 
the second metric is approximately flat). On the other hand, the first metric in the region 
of 5 ^ 1 experiences a power-law accelerated expansion, which is converted to deflation in 
the limit 5^0. Details of these solutions are given in Appendix C. 

Let us now describe the behaviour of our solutions obtained when working with the 
"nicer" potentials, able to confine the motion of 7, such as the modified Pauli-Fierz potential 
(|96|). In that case we numerically found that the presence of terms in the potential ^(7) 
quadratic and quartic in 7 are sufficient to bound the motion of 7 in a finite range. We then 
insert the solution for 7 of = into the evolution equations for the other "timelike" 

degrees of freedom, i.e. for a and (3, or equivalently for a = a + (3 and 5 = (3 — a. A 
numerical study of the evolution of 7, o and 5 then shows that the expected behaviour 
of the two vertically separated "particles" connected by a "spring" essentially holds true, 
with, however, an important effect linked to the "center of mass" motion. If, for instance, 
we consider the case where initially some "particle", say a M , is "below" the other "particle" 



one with 5 > 0, 7 < -2 - \ 



|5 and the second with 5 < 0, 7 > 2 — |5. Note that both 



Kln « 1). 
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we found (for the initial velocities a, j3 we considered) that the "lagging" "particle" a- 



starts to chase the "leading" one, and tries to overtake it. In other words, if 5 — (3 — a — 
i^niP^ — o^) is initially negative^, it will tend to increase. 

However, we found that the further evolution of 5 disagrees with the naive expectation 
of a timelike harmonic oscillator. Indeed, we found that 8 initially increases, but then sta- 
bilizes, after some damped oscillations, around a negative final value without ever crossing 
zero. [Correspondingly 7 also stabilizes, after some damped oscillations, around some final 
value.] This behaviour is clearly due to the coupling between the evolution of 5 and the 
evolution of a, i.e. the coupling between the "relative motion" of the pair of "particles", 
with its "center of mass" motion. The latter tending to damp the former, and to maintain it 
away from zero, i. e. away from a zero-separation configuration. We expect such a behaviour 
to be generic, though its details might depend on the chosen initial data. [For instance, 
for some well chosen initial velocities the lagging particle should be able to overtake the 
leading one so that one could have one or a few vertical oscillations before the particles 
stabilize into a permanently "chasing" configuration with non-zero vertical separation.] 

What is physically interesting in this behaviour is that we have here a locking mechanism 
by which the coupled evolution of isotropic metrics — e 2a//3 <5^, xjL — e 2 ' 3 / 3 ^, i.e. a 
"bi-Friedmann" configuration, locks itself up, after some damped oscillations, in a stable 
"chasing" configuration where 5 = (3 — a is constant, as well as 7. In view of the zero energy 

constraint (0), taken in the limit where <5 M = and o^o^ = — |cx 2 , the final configuration 
must be such that: 



Therefore, V must end up being positive, and a 2 ends up being a constant. In other 
words, the locking of the two metrics fixes the value of the bimetric potential to a positive 
value which then behaves as a usual cosmological constant term driving an exponential 



metrics. Note, however, that the physical values of the two expansion rates, as measured in 
each corresponding metric, are different: Hi = e 7//2 ifo and Hi = e~ 7 / 2 f/o where Hq = a/6. 
In view of the potential physical importance of such a "bi-de-Sitter" locked configuration, 
let us consider the general conditions for such a stationary configuration to exist. 

8 Remember from Fig.|] that 5 is opposite to the vertical separation between /3 M and a M because of the 
timclike character of n M . 




(102) 



inflation, at the same rate (as measured in the average proper time) a = (3 = |, for the two 



27 



Starting, for instance, from the full set of "timelike" equations given in Appendix A, one 
easily finds that the necessary and sufficient conditions for a locking solution (i.e. 7 =const., 

5 = (3 — a =const., a = = ~ =const.) to exist are: 

1 dVh,S) 1 , fj-5\ . . 

— tanh - (104) 



V 9 7 2 V 2 

with the additional requirement that V is positive so that (|102|) can hold. Then the constant 
common rate of expansion a = (5 = | is given by Eq. (|102|) . Note that Eqs.( |103|) , ( |104| ) 



represent two equations for two unknowns. One generically expect this system to admit 
solutions. If we consider the class of potentials that depend only on <j\ and a 2 , V = 
V(cri,a 2 ), using (|57j), (^8|) with r = 0, i.e.: 

<7! = 2( 7 + <y) , a 2 = 4^7 2 + ^ (105) 

one easily sees that the condition ( |103| ) simplifies to: 

There are two types of solutions of the above equation: either 5 = 37 or d a2 V = 0. 
Note that we already encountered the condition d a2 V = in Appendix B as the necessary 
condition for nontrivial "symmetry breaking" effects to occur (i.e. that V, considered 

as a function of the relative shift vector b, be roughly of the symmetry breaking form 
V = Vq — b 2 + b A ). However, the potentials that we are currently considering (including 
their A-modified versions, e.g. (p6|)) do not admit "critical points" where d a2 V = 0, if 
A > 0. Let us then concentrate on the other type of universal solution of the constraint 
namely: 

5 = 37 (107) 



Inserting this general solution in (104), we then get one constraint for one unknown 



(say 7). Note, however, that not all solutions of this constraint lead to locking solutions. 
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Indeed, one must further require the corresponding value of V to be positive so that one 
indeed ends up with a (real) "bi-de-Sitter" solution with common expansion rate: 

a 2 = $ 2 = - (108) 

4 cosh (2=S) V ; 

The solutions satisfying (|107|) lead to the following metrics: 



ds 2 2 = e +7 ds 2 



where ds 2 = -dt 2 + e a,z 5 tlv dx^dx v (109) 



Note again that the physical expansion rates are H± = e 7//2 ifo and Hi = e~ 7//2 i?o where 
Hq = a/6. One can check that such bi-de-Sitter solutions in bigravity admit, for general 
classes of couplings, generalizations to multi-de-Sitter solutions in multigravity theories 
with N coupled metrics. The conditions for this to happen is that = Aj/3, where 
Tj — Ti+i — li an d Aj = a i+ i — ctj, with i — 1, . . . , N — 1 (%, ati being the logarithmic time 
shifts and scale factors respectively for the i-th metric) and additionally that there exist 
real solutions to the N — 1 equations (similar to ( |104| )) that determine the IV Then one 
indeed finds that in the appropriately defined averaged proper time coordinate: 

ds 2 = e f 'ds 2 , where ds 2 = -dt 2 + e a/3 5^ u dx> l dx' / (110) 

i N 

with a = a i + and /j = Tk — Tk ■ 

k=l k=i+l 

For instance, for the original Pauli-Fierz potential, Eq. (|104 ) with ( |107|) leads to the 



constraint: 

7tanh7 = l (111) 

which has the two solutions 7 C ~ ±1.20. However, the corresponding value of V (V — 
— 2m 2 7 2 ) is negative. By contrast, we find that the modified Pauli-Fierz potential (|96|) 
with the quartic term \a\ admits locking solutions. In the latter case, Eq. ( |104j ) with ( |107| ) 
leads to the constraint: 

7tanh 7 = | A7 / 2 _ 1 (112) 

which has for any A > two pairs of roots (each pair consisting of two opposite roots), 
one pair with §A7 2 — 1 < (which tends to -y c ~ ±1.20 as A — > 0) and another pair with 
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|A7c — 1 > (which behaves as 7 C ~ ± yA= + |J as A — > + ). For A < there is only one 
pair of opposite roots. The value of the potential, on the other hand, is: 

V = 2mW c Qa 7c 2 - l) (113) 

which is positive only if A > and additionally \^1 2 — 1 > 0. Thus, the requirement of 
positivity of the potential excludes the A < case and keeps only one pair of solutions for 
the A > case. Note that any positive value A would give a locking solution. 

Summarizing, general classes of potentials admit timelike separated configurations which 
lock in a stable^ inflationary regime corresponding to a "bi-de-Sitter" configuration, with 
the same expansion rate. We see that bigravity can therefore play the role of dark energy in 
driving cosmic acceleration. But, contrary to the spacelike separated case (which led to an 
anisotropic type of slow-roll inflation, followed by power law expansion driven by spacelike 
oscillators) the pure timelike separated case is consistent with an isotropic de-Sitter-like 
inflation. When comparing such a final behaviour, in the modified Pauli-Fierz case, with 
the discussion of the original Pauli-Fierz case in Appendix C, it is striking to note how a 
simple (and natural) modification of the potential can drastically affect the set of bigravity 
solutions. This is a clear illustration of the concept of "universality classes" of bigravity 
potentials discussed in [p4fl . 



It is interesting to compare our locked solutions to the accelerating solution found in 



25. 26] within the context of brane- induced gravity. In the late time limit the solution 



of pq,p6[ is given by: 

ds 2 = dy 2 + (1 + \y\)dsl , where ds 2 = -dt 2 + e 2m b ^dx^ dx v (114) 

The fact that the four- dimensional part of the above five-dimensional metric is conformal 
(with a y-dependent conformal factor) to a fixed de-Sitter metric is similar to the particular 
case ( |110| ) of our general locked solutions. Eq. ( |114j ) can be easily understood as following. 



The five-dimensional Einstein action contains the term J \ / —g[tr(g~ 1 d y g) 2 — (trg~ l d y g) 2 }. 

This term can be viewed as the continuum limit of an infinite sum of "nearest neighbour" 

interactions Yl^{9i+i9i) ( see till) wnen considering the N — ► oo limit of multigravity 

(where N is the number of coupled metrics). Based on this reinterpretation of the particular 

9 We have not investigated the general conditions for stability of such locked configurations when they 
exist. However, they seem to be numerically stable even when introducing a tilt angle away from the 
vertical direction. 
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five-dimensional model studied in [^,[26|], we can view it as the N — > oo limiting case of 
the general class of de-Sitter locked solutions mentioned above. 



6 The brane motivated potential 



The brane motivated potential flT8|) with the previous notation reads: 



V(j, 5, r) = m 2 



C ° Sh ' 2vlV ^ 2 + 3~^ 2 + ^ ~ 



cosh 



5 + 7 



;ii5) 



6.1 Absence of timelike separated worldline solutions 

If we look for solutions where the two worldlines are purely timelike separated (r = 0), the 
Lagrangian (B3) becomes: 



£ = e^ 2 



■i cosh CL-f) (a 2 + 5 2 ) + I sinh {^A ) &S _ y ( 7j *) 



(116) 



with the potential: 



V(-f,S) = m 2 



cosh I — — —I — cosh ^ 



(117) 



In this case, one finds that the two constraints (0), (|49|) force 5 = 7 = 0. Thus, the 
only solution is & = 0. In other words, there are no nontrivial solutions with r = 0. Note 
that this does not preclude that the existence of timelike separated solutions with non- 
vertical separation, i.e. with 5 M timelike but not parallel to n M . Note also that there is no 
contradiction between the facts that: (i) the brane potential reduces to the Pauli-Fierz one 
for small 7 and (ii) the Pauli-Fierz potential admits purely vertical timelike evolutions, 
and (iii) the brane potential admits no such evolutions. Indeed, the Pauli-Fierz timelike 
evolutions only exist when 7 is of order unity. 



6.2 Spacelike separated worldline solutions 

In the opposite regime where the two worldlines are required to be purely spacelike separated 
(5 = 0) one again finds that the only way to freeze 5 to zero is to have also 7 = 0. Then 
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the Lagrangian ( |65|) becomes: 



C = Qr 2 - h 2 - V(r)j 



;ii8) 



with the potential: 



V(r) = m 2 



cosh I — r ] — I 
3 



(119) 



Let us note here that the exponentially growing structure of the above potential (with 
a coefficient of r of order unity), makes us expect that there will be no slow-roll region. 
We can (at best) expect that such an exponential "spring" will lead to power-law inflation. 
Asymptotically, for r ^> 1 the action can be written as: 



S 



1 . a m s/3, 

2 a z ' 

3 3 2 



dte a/2 I -/■' - -cT 



1 

— ( 

3 



where we have changed our time variable to dr 



dT ( 2 r '2 _ 4^2 



dt 



(120) 



^72. Defining new variables as: 
= — a + r 



(121) 



we can have a considerable simplification in our action: 



[122) 



This "Toda model" is integrable because the £2 variable is ignorable. We have: 

6 = 6o+2>r (123) 

where p is the conserved momentum of £2 and £20 an integration constant. From the 
Hamiltonian constraint we immediately obtain for p 7^ 0: 



£1 = - log 



7m 2 . 2 /pr 
— — smh — 
6p 2 V 2 



(124) 



where we have set an integration constant, which accounts for time shifts, to zero. For the 
special case of p = we get: 

~7m 2 T 2 ~ 



f 1 = - lof 



24 



(125) 
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There is an obvious symmetry in the simultaneous reflections r — > — r and p — > —p. 
These map solutions of expanding volume (increasing <r) to ones of shrinking one. We will 
be interested in expanding volume solutions where r < 0. The parameter p represents a kind 
of initial kinetic energy. Let us focus here on the simplest case p = 0, which corresponds 
to the usual rolling-type inflationary solution with kinetic energy small compared to the 
potential energy. See Appendix D for a discussion of the general case p 7^ 0. 

For this case (p = 0) we have for any value of the time parameter r: 



a 



18 

7 



log \ mr \ + const. 



V2 



-a + const. 



(126) 



The proper time t (which is proportional to the individual proper times t%, £2) is related to 
r via t oc r~ 2 / 7 , so that the solution ( |126| ) reads in proper time: 



a = 91og(mi) + const. , r = — 3v^2 log(mt) + const. 



(127) 



By rescaling coordinates to absorb the above constants we finally see that, in conformity 
with naive expectations, each metric has a power-law behaviour, and if we parametrise it 
in the standard way as: 

3 



ds 2 = -dt 2 + J2( mt ) 2p "( dx 



we have that the Kasner exponents p M are: 



:i28i 



+ 



and 



P(2) 



+ 



(129) 



We should stress here that, when we bring one of the metrics to the above form (|128[) by 
absorbing the various constants in redefinitions, the second one will not have the form 
128|) . Thus the exponents we have written above are in the coordinate system where each of 



the metric (but not both) simplifies. The above exponents do not satisfy the quadratic (zero- 
mass-shell) Kasner relation, since X^G° M ) 2 — (j^^P 1 ^ = ^ n addition, J^^P^ = |; 
which means that each metric's volume Vi expands as a function of its respective proper 
time as V{ oc , with % = 1, 2. Note that the latter volume expansion has an accelerating 



33 



behaviour, but as expected, much less pronounced than the case of the polynomial Pauli- 
Fierz potential which led to exponentially inflating solutions for such a rolling initial state. 
The evolution is highly anisotropic, with the anisotropy parameter being: 

A x = A 2 = 1 (130) 

In all the above discussion it is always assumed that r ^> 1 when the simplified action is 
valid. For rCl the potential is identical with the Pauli-Fierz one and the r, o will behave 
exactly as we have discussed in the previous section for this case. 

Let us now check again if the obtained solution is stable against perturbations of S. We 
have to recalculate the quantities A and B of Eq.(^) and their limits for r ^> 1. They 
read: 

A = --e^ 2 , B^ Z ^—e°l 2 (131) 
33 297 y ' 

Thus, the variation of the action simplifies to: 

5S = [ dte°' 2 (--6 2 - ^-5 2 ) (132) 
J V 33 297 J K J 

From the extremization of this action, we get the following motion for 5: 

(Vt ,. s I -\!i-mH 2 \ - C,K; 



1 





(133) 



where I n (x) and K n (x) are the modified Bessel functions of first and second kind respec- 
tively. Hence the solution is unstable at late times. This instability is again linked to a run 
away of 7 towards large values (for which the above perturbative treatment breaks down). 
Numerically we observe that as soon as I7I starts to increase, it tends to run away very 
quickly towards infinity. The behaviour in this limit is similar to the one presented for the 
Pauli-Fierz potential. 

Having understood that the brane potential ( |115| ) is, like the Pauli-Fierz potential, 

not efficient enough^ in confining 7 in the long term, we have also numerically studied 

10 This is a priori a bit surprising because Vt, rane (7) increases exponentially as I7I — > 00. However, 
VbraneM increases no faster than the kinetic terms in Eq.(B5|), i.e. like ae 1 ! 2 + 6e~ 7 / 2 , and the asymptotic 



coefficients a and b depend on the value of S. Therefore, it it can become subdominant with respect to the 
(sign-changing) kinetic term contributions. 
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the cosmological evolutions obtained with a brane potential modified by the addition of a 
term oc Acrf . Our (partial) numerical study then shows that such a better-confining brane- 
type potential leads to long-surviving solutions which are qualitatively similar to the ones 
discussed above (in the modified Pauli-Fierz case). Again, we find that the simple picture 
of two "particles" connected by a "spring" is a good qualitative guide. We found spacelike 
separated solutions which start, as above with a power- law anisotropic expansion, and which 
end up in some (power-law expanding) "spatially" oscillating regime. The modification also 
allows now for purely timelike separated solutions which end up, after "locking" , in a bi- 
de-Sitter final configuration. 



7 Transition between matter domination and vacuum 
bigravity domination 

In the previous sections we have seen that, for both potentials that we have studied, there 
exists a period of acceleration for one or both of the metrics. It is tempting to suggest that 
this acceleration due to the metric coupling potential can be identified with the current 
acceleration of the universe. Let us see if this type of "dark energy" can be phenomenolog- 
ically acceptable. A basic ingredient of any theoretical model of dark energy is that there 
exist a smooth transition between a matter dominated universe and a dark energy domi- 
nated one. To discuss this transition for our bigravity system we consider the addition of 
matter densities in the action (p8|), having a scaling law p\ m = pi m oe~ a for the first metric 

and p2m = P2mo£~ 13 for the second one, where pi m o and p2mo are some initial values of the 
two energy densities: 

C = e a ~^ a»a v + e^ 2 - e^+^+^V^, - e> lm0 - e> 2m0 (134) 
The two Friedmann equations (in the 7 = gauge) are modified as: 

e f af 1 a (t + p x + p lm = (135) 

e" 7 /?^ + P 2 + P2m = (136) 

where, keeping the same notation as in the previous sections, we denote by p\ and pi the 
massive graviton dark energy densities. 

We have done some numerical experiments on this transition, notably in the case of the 



quartically- modified Pauli-Fierz potential (pi). In the examples we studied the confining 
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power of V(7) was strong enough to keep 7 in some finite range, even when the potential 
contribution in Eq.( p.34j ) is very small compared to the other terms. Then we have observed 
numerically that there was a smooth transition between matter domination and dark energy 
domination for both the spacelike and the timelike worldline separation cases. 

Rather than describing in detail these numerical results (which are not particularly 
illuminating) let us explain with the help of simple analytical formulas what, we think, 
is the general simple behaviour of this transition in bigravity theories. Note first that 
several types of transitions can occur. If we are interested in the phenomenology on the 
"first brane", g\jj^, either pi overtakes p\ m before p<i overtakes p2 m , or the reverse. Our 
(partial) numerical experiments indicate that these different transitions lead to the same 
phenomenological features on the first brane. For simplicity, let us then consider the case 
where p\ overtakes p\ m before p2 overtakes p2 m . Then, during matter domination (on 
both branes), both an d x\v expand as t 2//3 (as 7 is fixed the average proper time 
t is proportional to both "intrinsic" proper times). This means that x]i V {t) = a fivt 2 \ 
xfiu(t) = b^t 2 ' 3 . Here, a M „ and b^ u are some constant tensors. Each one can be reduced 
(by a change of spatial coordinates) to the normal (flat) isotropic Friedmann form, say 
a' = 5^. However, it is important to realize that the spatial transformation = A^x ,l/ 
which will "isotropize" x\ v i wm , * n general, leave x\ v m a general, apparently anisotropic 
form V 7^ XS^. In fact the "difference" between and b^, i.e. more precisely the matrix 
(Xi 1 X2)u = aiLK ^K,v can a priori take arbitrary values. The corresponding logarithmic 
eigenvalues pa = log(eigenvalues(a _1 b)) essentially define an "initial" field-space vector 
5^ = (3^ — a M separating the two "particles" during matter domination. In other words, 
in the mechanical language of the model of Fig.|l], during matter domination (where the 
"spring" between the two worldlines has only a subdominant effect^]) the separation vector 
between the two "particles" remains constant: <5 M = /3 M — a M =const. Under our assumption 
about the ability of ^(7) to confine 7, we therefore see that, during matter domination (or 
for that purpose, also radiation domination, where X/LW = a ^ 1 ^ 2 ^ e ^ c the effective 
vacuum energy seen by each metric, e.g. p\ = e ( ^ 7+ ^ //2 (iV A (7, 5^) — 7 V(7, S^)), remains 
essentially constant. It plays therefore the same role (modulo possible 7-modulations) as a 
cosmological constant, during any period where V is subdominant. 

11 Note, however, that, even when V is negligible with respect to the other terms, its presence remains 
crucial for pinning 7 to some specific value, which we assume here to remain 0(1). 
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However, when, because of the decrease of p\ m oc e~ a during expansion, p\ starts dom- 
inating over pi m the situation will change and the separation vector 8^ will start to vary. 
In fact, if we are in the case where the "initial" value of p\ (as determined by the nearly 
constant values of 7 and 8^ during matter domination) is positive we simply expect that p\ 
will take over p\ m the role of driving the expansion. Such a situation can a priori occur for 
any "spacetime direction" of the connecting vector 8^. It is here that a bigravity origin of 
dark energy ("tensor quintessence") can lead to interesting phenomenological predictions 
that differ from the standard (scalar-quintessence-like) models. 

Two contrasting cases can occur, depending on the timelike versus spacelike character 
of the initial connecting vector 8^. First, let us consider the case where 8^ is timelike, and 
even, for simplicity, of the extreme timelike type 8^ oc Sn^, i.e. r = 0. In this case, if 
there exists a stable locking configuration (as described above), we expect 8 to ultimately 
stabilize (after some oscillations and after the second brane transits to a vacuum dominated 
state) in this locked state: 8 = 37 =const. Note that, in the intermediate stage where the 
second brane is still matter dominated there will be an "external force" acting on 8 = (3 — a 
(through the equation of motion of (3) so that one expects 8 to deviate, in a time- dependent 
way, from the final locked configuration. In other words, we expect in this case the effective 
equation of state of dark energy: 

p i V - 2d s V f 

Wl = jr -V-iKy (137) 

to depend on time, before ultimately stabilizing into the standard "vacuum" value w\ — — 1 
corresponding to the final, bi-de-Sitter locked configuration. 

Let us now consider the other possible case where 8^ is spacelike, and even, for simplicity, 
of the extreme spacelike case where 8^ _L n 11 , i.e. 8 = 0. In this case, an interesting 
prediction of bigravity is that the expansion rates of each universe (which, during matter 
domination were both intrinsically isotropic) will start to become anisotropic. Indeed, we 
expect that (at least after the second universe also transits to a vacuum dominated state) 
the solution for the connecting vector 8^ (which started with the initial conditions 8^ = 8$ , 

8^ = 0) will quickly be attracted towards the spacelike slow-roll solutions described in 
Section 5. This means that, within one Hubble time or so after the transition, the intrinsic 
expansion rates along the three special spatial directions (corresponding to the eigenvectors 
of the matrix a _1 b) will differ from each other. We expect, however, like in Eqs.(jS^)- (|85|) 
above (which are not directly applicable because they were derived for the original Pauli- 
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Fierz potential) that, if the magnitude r = y|A^A M of the (purely spacelike) connecting 

vector <5 M = A^ is large compared to one, the spatial anisotropies in the rates of expansion 
will be rather small. As for the time dependence of the expansion, it will be initially of the 
slow-roll inflationary type (as in Eq. Q7^),(|5Q)), but will ultimately end up as a power-law- 
type expansion as r becomes 0(1) and starts oscillating around zero. 

We have checked the occurrence of the two types (timelike versus spacelike) of transition 
to bigravity vacuum dominance in some numerical experiments. We have, however, not tried 
to make exhaustive experiments because there are many parameters that can influence the 
actual time- dependence of the transition: the shape of the potential, its ability to confine 

7 or not, the value of the the ratio between p lm and p 2m , the presence or absence of a 
"locked" solution, etc. It is too early to embark on such an exhaustive study. Our main 
purpose in this work was only to delineate some of the qualitative new features of possible 
bigravity origin of cosmic acceleration. 

8 Conclusions 

We have explored general spatially-flat cosmological solutions (of the anisotropic Bianchi 
I type) of classes of nonlinear bigravity theories. Even within this restricted class of ho- 
mogeneous cosmologies we focused on special cases. We did not explore the possibility 
(mentioned in Section 3) where, due to a type of spontaneous symmetry-breaking mecha- 
nism, the relative shift vector 6 M = e~^(&2 — &i) be non-zero. After a brief discussion of the 
structure of the evolution system for two coupled spatial metrics, we restricted ourselves 
to the simple case where the two metric tensors can be simultaneously diagonalized. Even 
this simplified case leads to very rich dynamics which can be conveniently described in 
terms of a mechanical model (represented in Fig.|I]): two "relativistic particles", moving in 
a (2 + l)-dimensional Lorentzian space, and connected by some non-linear "spring", i.e. 
interacting via some bigravity potential V12 = (fi , ifi , 2) 1 ^ 4 ^(gjf 1 g2). 

One of the first important conclusion of our study is that the long-term behaviour of this 
coupled system crucially depends on the ability of the potential V to confine the evolution 
of the relative lapse 7 = log (N2/N1). Due to the former "gauge nature" of 7 (when the 
potential is absent), i.e. the absence of kinetic terms for 7, the equation of motion of 7 
is algebraic. We found that the continued existence (in the long term) of a solution for 
7 sensitively depends on the nature of the function ^(7). For instance, we found that 
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the potential V(gi 1 g2) derived from the five- dimensional brane constructions has only a 
marginal ability to confine the evolution of 7 to a limited range of variation. Indeed, we 
found that many (and maybe most) solutions of bi-cosmology, with such a brane potential, 
evolve, after some finite time, into a state where 7 quickly runs away towards infinity. 
By exploring the behaviour of physical observables near the moment where I7I — ► 00, we 
have shown that the run away does not correspond to any observable singularity in either 
of the two metrics. We have argued that this run away only signals a breakdown of the 
effective four-dimensional description that we use. Indeed, it seems that, as I7I — > 00, some 
previously heavy modes become light and should now be taken into account in the effective 
action. 

Having understood the root of this run away behaviour, we have focused our physical 
study of bigravity on the class of potentials ^(7) which have the strong-enough confining 
property with respect to 7. A simple example of a potential in this class is the "quadratic 
plus quartic" Pauli-Fierz-type potential (j2ip . Such a potential allows for solutions which 
evolve on long-time scales, without encountering any breakdown linked to a run away of 
7. [Note, however, that we do not claim that this is true for all solutions. It is certainly 
possible to concoct initial data leading to a 7 runaway after a finite time.] 

When using such 7-confining potentials (or when considering, as we do in most of the 
text, the effect of any potential up to times smaller than the moment of quick 7 runaway) 
we have found that the qualitative behaviour of generic cosmological solutions can be nicely 
understood in terms of the mechanical model of Fig.|l| For instance, when the separating 
vector 5^ = f3^ — a^ between the two "particles" is spacelike, and the potential is polynomial 
in <5 M (and attractive, as the modified Pauli-Fierz potential (|21|)), the coupled motion of 
and /3 M is similar to slow-roll inflation. The separating vector <5 M plays the role of the inflaton, 
and drives an exponential-type expansion of the vertical position of the "center of mass" 
|<7 = + a) (which represents the average volume of the two metrics). A qualitatively 
new feature of this type of bigravity slow-roll inflation is its growing anisotropy. As the 
connecting vector gets smaller, each metric tends to expand more and more differently 
in three spatial directions (linked to the "direction" of the vector 8^). This anisotropic 
slow-roll inflation ends up in a regime where the "spring" connecting the two "particles" 
makes them oscillate along a spacelike direction (see Fig.|3|). Similarly to the oscillatory 
period following slow-roll (for a chaotic inflation type potential, e.g. V((p) = |m 2 2 ), 
these bigravity oscillations lead to a power-law expansion law. We expect such anisotropic 
accelerating solutions to exist in a general multigravity theory, and thus also in the brane- 



39 



induced model fl9, 20|, which, as we said above, can be viewed as the N —>■ oo limit of the 
particular "nearest neighbour" interactions multigravity model. 

When the separating vector 5^ = (3^ — is timelike, we found (when using the 7- 
confming potential fl2T|)) a remarkable phenomenon of "locking" of the two metrics. [We 
have indicated in Eqs.( |103| ), ( |104| ), ( |108| ) the general conditions under which this "locking" 
phenomenon occurs.] In the visual language of Fig.|I], the two "particles" lock in a perpet- 
ual "chasing" configuration where their vertical separation tends to a non-zero constant, 
while their "center of mass" continues to move upwards. In bimetric language, the locking 
corresponds to a bi-de-Sitter configuration: each metric expands exponentially, and the two 
expansion rates are equal (in the average proper time). Contrary to the spacelike case where 
anisotropies played an important role, here this configuration is obtained for isotropically 
expanding metrics. It would be interesting to explore the basin of attraction of such a 
locked state among generic (timelike-type) bigravity evolutions. As indicated above, these 
locked configurations admit (provided some system of N — 1 equations for N — 1 unknowns 
admit a real solution) a multi-de-Sitter generalization in a general multigravity model. [As 
we mentioned, this is a way of interpreting the solution of [25|, 261-1 Since we found that the 
bigravity locked solution was stable, we expect this feature to extend to the multigravity 
case. 

From the phenomenological point of view, one of the major conclusions of this work 
is that bigravity cosmologies generically exhibit a period of cosmic acceleration for one or 
both of the metrics. This conclusion applies even to the case of "bad" potentials which 
cannot permanently constrain the evolution of 7 to a bounded range. This result suggests 
that bigravity could be the origin of the observed cosmic acceleration, i.e. that it could be 
the the source of dark energy. In other words bigravity naturally defines a kind of "tensor 
quintessence" . In brane models, the mass parameter m 2 appearing in the potential V is an 
exponentially decreasing function of the interbrane distance. It is therefore not unnatural 
to have an m 2 as small as it is required to explain the observed cosmic acceleration (i.e. 
m ~ 10" 33 eV)Q. 

Our preliminary studies of the transition between matter domination and vacuum dom- 
ination seem to indicate that (at least for classes of potentials) this transition can be as 
smooth as in the usually considered dark-energy models (such as a cosmological constant, 
or some type of scalar quintessence). It is, however, interesting to note that, at least in 
12 Let us note that the parameter fi appearing in (0) would then be [i ~ 10 _3 eV. 
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the spacelike separated case, bigravity makes qualitatively new predictions: it predicts a 
growing anisotropy of the expansion of the universe. It would be interesting to study the 
imprint of this phenomenon (which started to take place only "recently" , i. e. for redshifts 
z < 0.5) on observable phenomena, and notably on the Cosmic Microwave Background. 

Finally, on a more speculative view, it would be interesting to explore the possibility 
that bigravity explains the primordial inflation needed to explain the gross features of our 
universe. For this, one would probably need a mass scale of order m ~ 10~ 6 Mpi. If we 
contemplate a "spacelike" scenario, the needed large initial value of <5 M to have a long stage of 
bigravity slow roll inflation, might be naturally provided by the recently discovered generic 



chaotic behaviour taking place in (bulk) string/M cosmology [27,28,^9[. Indeed, the chaotic 
behaviour naturally leads, near t ~ t s trmg to very large "oscillations" in the (logarithmic) 
scale factors of the metric (considered at some spatial points). When comparing two 
metrics (either at two different bulk points, or on two branes) it is then natural to reach 
large values of (3^ — oA On the other hand, if we contemplate a "timelike" scenario, the 
bi-de-Sitter locked configuration might naturally explain primordial inflation. In this case, 
one still needs an exit mechanism (which could be provided by some instability linked to 7 
in the case where the potential ^(7) cannot indefinitely succeed in confining 7). 
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Appendix A: Equations of motion and constraints 



The Hamiltonian of the system is: 



n = e^ 2 



3 f *y $ \ ■ -) - } 

- cosh I — ^— ] (-p CT - p s + p s+ + + p A+ + p A Je 



+^ sinh (-^y~ ) 1 ~ ' 1 - i + )f ' + 1 



(A.l) 



where the canonical momenta are defined as: 



PA ± = r 



a/2 



v. - Y 



cosh ( — - — I A± — sinh I — - — I E± 



cr/2 



cosh I - - 1 5 — sinh ( - - I a 



cosh I ^— — ] a — sinh \ ^— — ) r) 



(A.2) 
(A.3) 
(A.4) 



The equations of motion in the "averaged proper time" t for the full Routhian are: 



• 5 equation 



2 „io d 

g <t/2 

3 eft 



e°"/ 2 ( cosh ( ^— — \ 5 — sinh ^ 



dv dv 



• (T equation 



-a/2_ 



dt 



e a/2 ( cosh 



a — sinh | ^— — ) 5 



(A.6) 



A± equations 



3 -<r/2_ 



ps± tanh 



7 



2 J 3 cosh 



-A, 



dV 



OA 



± 



(A.7) 
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3(p! + +p!_) _,tanh(^) 



-a/2 



+ 



cosh (^) 2cosh 2 (^)^ + 



(p s+ A + + Ps _A_) 



+ 6 Smh (^~ 



+ A 2 
cosh 2 pHj 



+ a 2 + ^ 



1 /'y-6\ ■ dV 

-3 cosh (V)^ + a7 = (A ' 8) 



In the above we have used the 7 constraint to simplify the 5 equation of motion. We 
additionally have the Hamiltonian constraint H = 0: 



3(p| + +p|_)e- CT 1 

+ - cosh 



/ K r — CUSll 

4 cosh (2=«) 3 V 2 



' A 2 , + A 2 

.cosh 2 pH) 



- a 2 - 5 2 



+ jjsinh (^y^ ) <t<) V = 



(A.9) 



The equations of motion in the "averaged proper time" t for the case where ps 4 
= are: 

• 5 equation 



dt 



e^ 2 ( cosh [ 5 - sinh ( ^ , a 



dV dV . . . 

97 + w (A ' 10) 



a equation 



,-<r/2_ 



e*/ 2 ( cosh I 2^1) b - S i n h fc* U 



V 



(A.11) 



• r equation 



3 rftlcosh(^) 



3p2 ff , fi-S\ dV 

-- *-f e _ff cosh -i + — 

2r 3 V 2 / dr 



s- ( A - 12 ) 
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• 7 constraint 



8r 2 



+ a 2 + 5 2 



cosh 2 (iz*) 



]cosh(l^)^ + ^ = (A.13) 



In the above we have used the 7 constraint to simplify the 5 equation of motion. We 
additionally have the Hamiltonian constraint H = 0: 



45 e ^ cosh (V) + r osh ( 1 2 i 



• 2 r2 

a — 



cosh 2 (^) 



+ ^sinh ( ?--A ) &5 + V = (A.14) 



Finally, the equations of motion in the "averaged proper time" t for the case of r = 0, 
when we have isotropic metrics, are: 



4 d 

3 It 

4 d 
3 eft 



e + 2a 



e p 2 







e 2 ( y 

c?7 85 

1 , . 8V 8V 

= e~ V H 1 

87 85 



(A.15) 
(A.16) 



We additionally have the two Friedmann constraints (which we have also used to simplify 
the above equations of motion): 



-a 

3 



_i (V 8V 
(V 8V\ 



(A.17) 
(A.18) 



44 



Appendix B: Perturbation analysis around the = 
shift vector in the cosmological metric ansatz 



In this Appendix we study under what conditions there might exist more solutions than 
the "trivial" solution 6 M = discussed in the text. The eigenvalue problem: 

(gftiM = (B.i) 

considered for 6 M 7^ but small, can be related (to 0(b 2 )) by standard techniques to the 
unperturbed one (denoted by an overbar) as: 

K = A a fl + M£p) (no sum) (B.3) 
where b a = e^6 M . The variation of the action when making a perturbation is: 

I Aa 



^El^V^o -^/)^ (B.5) 

~^ — A a 



and hence: 



^-^^ (r(2) ° - T(2)/) (n ° SUm) (R6) 



As said in the text, the equation of motion of b a is the vanishing of flB.6|) . There are only 



two possible ways to make the above quantity vanish: either b a = or T^ 2 \° — T a a = 0. 
The first possibility means that b a = is an isolated solution. Therefore, it is the second 
possibility which signals the threshold for the existence of new solutions, besides the trivial 
one. [We are assuming here that, as some parameters vary, the nonperturbative solutions 
can be made to coincide with the perturbative one.] On the other hand according to 
©,(0): 

T ( 2 ) o O _ T{2) a = _ 2e -^ {d ^ v _ d ^ v) = _4e-^/ 4 ( yUo - fr)d n V (B.7) 
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for potentials of the class V = V(cri,a2). Therefore, a general necessary condition for the 
possible existence of non-perturbative solutions is that d a , 2 V = admits solutions. 



Appendix C: Analytic solutions for the Pauli-Fierz po- 
tential in the timelike worldline separation limit 

In this Appendix we sill study an analytic description of the solutions of the extreme 
timelike worldline separation (i.e. r = 0) for the original Pauli-Fierz potential (|TB|). As 
discussed in the text, the solutions exhibit an initial stage of acceleration for large 5 and 
a period of deflation as 5 — > for the first metric, while the second metric remains ap- 
proximately flat. We will split up the analysis of this system into the two above-mentioned 
asymptotic regions. 

Let us note the equations of state for the two metrics obtained by (|70D, ([73]): 

"» = £ = - 1 + 2 P + M + M 7 (a2) 

• The cosmological initial singularity limit where 7 — > — oo and 5 — > 00 

We have already qualitatively described in the text the behaviour of the system in 
this limit. Due to this motion in field space, we can see numerically that for all initial 
conditions, the first metric experiences accelerated expansion. As an example to illustrate 
the behaviour of the system in this limit, we will consider the case where 7 initially lies very 
near the 7 = —2 — 15 line. Then, the solution follows this line in a very good approximation 
both back in time towards the initial singularity, as well as forward in time until 5 ~ 0(1). 
From the (3 Friedman equation we get that (3 ~ and from the a constraint that: 

a 2 = -3m 2 ( 7 + 2)e-^ +3 ) (C.3) 

But approximately at = —5 = 37, so we have a differential equation for 7. Integrating this 
we get: 

Erf(V-(7 + 2)) = l-^- (C4) 
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7 



Erf^ 1 ) 1 - 



tr 



log 



(C.5) 



where t cr = ^r, Erf(x) = ^7= J* e y dy is the error function, Erf( the inverse error 

function. The limit we have indicated is at t —>■ 0. From the properties of the error function 
we have that if t — > + then 7 — > — oo. In the other limit that t — > t cr we get that 7 — > —2. 
However, we never reach the latter limit, since close to that, our approximation breaks 
down because 5 ~ 0(1). In the region that this approximation is valid, we have: 



6 = -3( 7 + 2) , a 



-5 + 2C 



a 



-5 + C , P = C 



(C.6) 



where C is an integration constant. We see that the first metric expands, but one should 
describe this expansion in this metric's proper time. Asymptotically, for t — > we have: 



dtt = e~ 1,2 dt tx = 2y/t^i 
On the other hand the proper time of the second metric is: 

2 



dU 



z*l 2 dt 



to 



3Vt 



f3 /a 



(C.7) 



(C.8) 



Thus, the first metric is intrinsically inflating with the scale factor behaving as e a ^ 3 ~ tf, 
while the second metric is approximately flat. 

The equation of state ( |C.1|) , in this limiting case which we are examining, is for the first 



metric: 



Wi 



2 
3 



(C.9) 



3 7 + 2 

while for the second one, (]C.2[) leads to w 2 — > —00, since the 7 = — 2 — |5 line is the root 



of the denominator. 

Let us now check if the obtained solution is stable against perturbations of r. The 
variation of the Routhian in quadratic order is: 



5C 



W2 



r 



3" \^ cosh (2=2) 3' 
From the extremization of this action we get the following motion for r: 



6 V 2tt 



(CIO) 



(C.11) 
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which is growing only logarithmically as t — ► 0. Thus, the solution is stable in very good 
approximation. This is in agreement with our numerical study. 

• The final state limit where 7 — > — 00, 5 — > 

The final state of the evolution of the system is independent of the initial conditions. 
As 5 — > we find that 7-^—00 and the first metric deflates and becomes asymptotically 
flat in infinite proper time, while the second experiences a finite (in proper time) period of 
inflation. The runaway of 7 in this limit is an unavoidable fact of the Pauli-Fierz potential, 
as we have discussed in the text, because the action ceases to have an extremum at finite 7, 
when one of the two worldline velocities tends to zero. In order to study this limit we need 
to do a different approximation to the equations of motion. Combining the two Friedman 
constraints (|100|) , ( |101| ) and keeping leading terms we get the following relation: 



The equations of motion on the other hand are approximated by: 



Ad_ 

3dt 




r 



7 a+£ 



(C.13) 




(C.14) 



Using (|C.12|) we have: 




(C.15) 



d_ 

dt 




(C.16) 



Then since a(3 — ► 0, we can neglect this term and obtain the system: 




(C.17) 



a — — 





(C.18) 



48 



From this system we get: 



' a 



(C.19) 



which putting back into the second equation gives: 



« = --log --1 



This can be solved in the region where a — > C\ and gives: 



(C.20) 



a 



l + li^ 



Ci 



m 2 . . , / m 2 , 

(t - to) log ( 7^r\t-t \ 



2d 



2Ci 



(C.21) 



where li(x) = is the logarithmic integral and li*- ^(x) its inverse function. The limits 

which have been used are that for x — > + , li 

This shows that a diverges to — oo. Asymptotically, the function a is: 



. ^ y x) — > and for x — > , li 1 - 1 - ) (x) — > x log |x| . 



a « C 2 + Ci(* - t ) - — (* - to) 2 log ( — |t - t | 



The asymptotics for the (3 function is: 

■m 2 

^ = C 2 + — (t -to) 2 log 



2 / 9 

m ,. . s2 / to 



2Ci 



* — *o 



and for <5: 



On the other hand 7 is: 



- .. ( — 1 

a 



5 = -C 1 (t-t ) 



7 



2 log 
2 log 



777-^ /" ?77, ^ 

■(*-*o)log H7rl*-*o| 



2d 



2d' 



(C.22) 



(C.23) 



(C.24) 



(C.25) 



The proper time in the first metric is: 



2Ci. 
— log 



Ioe 



to 
2d 



(C.26) 
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thus, the singularity point to hi the proper time of the first metric is at t\ 
other hand, the proper time for the second metric is: 



oo. On the 



dt 



2 = e 



^2 — ^20 + 



2d 



m 



2d 



{t-ta 



log 



m 



2d 



\t-U 



(C.27) 



and it is finite at t —>■ to- 

Thus, the first metric is intrinsically deflating with the scale factor behaving as: 



e ~ 



2 i 



1 — — -e~' "' 
3m 2 



(C.28) 



while the second is intrinsically inflating as: 



(C.29) 



The equation of state ( |C.1| ) on the first metric, in this limit, is: 



w 1 = 2- 





oo 



while on the second one (|C.2|) leads to the opposite effect: 



(C.30) 



w 2 = -2 



— oo 



(C.31) 



Appendix D: The general p ^ case evolution for the 
brane motivated potential 



In this Appendix we will discuss the evolution of the system for the brane motivated 
potential for the general case where the "initial kinetic energy" p is non-zero. The solution 
of the equations of motion is given by (|123|), (|124[). The initial "incoming" solutions are 



then for r < where the kinetic energy of the system is large in comparison with the 
potential energy Then we can write the solutions of a and r using (|121|) and separating 
the various constants (and subleading terms) as: 
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For p > 



a — > - (3 - V2)pr H , r — > a + 



(D.l) 



• For p < 



a ^ --{3 + V2)pr + ■ ■ ■ , r -> -<r + 



(D.2) 



On the other hand, the final state solutions are for r — > when, due to the expansion of 
the two metrics, the kinetic energy of the system has become subdominant in comparison 
with the potential energy. Then we have asymptotically for both signs of p: 



18 



log \pr\ + 



v/2 



(D.3) 



Note that the solutions in the latter epoch, have the same scaling law as the ones for p = 0. 
In the above language we see that in the case where p > 0, r initially increases until a 
maximum value and then shrinks to zero. On the other hand for p < r always increases. 
In all cases a increases and thus the volume of each metric expands. 

Now, we need to go back to proper time t to see the behaviour of our solutions. For 
r C both cases have a = cr + ■ • • with different constants c, whose value does not have 
any significance as we will see in the following. Then the proper time and the expression 
of a as a function of it are: 



t 



a = 2 \og(pt) + 



(DA) 



where again we ignored unimportant constants. By rescaling coordinates to absorb the 
latter constants we finally see that each metric has a power-law behaviour, and if we 
parametrise it in the standard way as: 



ds 2 = -alt 2 + ^{pt) 2p " \dx»f 



(D.5) 



we have that the Kasner exponents p M for p > are: 



1 + J - 

3 + V3 
1 1_ 
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and 
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1 + J - 

3 + V3 
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3 



(D.6) 



51 



and for p < they are the same with a flip on the sign of the second addendum of 
the first two exponents. The above exponents satisfy the usual quadratic Kasner relation 

X^G° M ) 2 — = 0> as well as YlfiP 11 = 1> which means that each metric's volume 

Vi expands as a function of its respective proper time as v oc t iy with i = 1,2. Thus, this 
evolution is highly anisotropic with: 

A 1 = A 2 = V2 (D.7) 



One can notice at this point a potential paradox because an exact Kasner metric is 



known to be "on the light cone", i.e. to have a^a^ 



W 2 -r 2 ) 



0. On the other hand, 



from the Hamiltonian constraint (|7q) the same quantity should be very large for r> 1. 
This can be understood if we include next to leading order terms in our asymptotic solution. 
These will modify the a asymptotic (|D.4[) by a term linear in t and also the r asymptotic 
with a term linear in t with different coefficient. This will immediately render a 2 — r 2 very 
large as expected. 

For the case where r — > 0, we can express the proper time and a as a function of it as: 



t = Ct-~ 



a = 9 log(pt) + • 



(D.8) 



Thus, the exponents this time are: 



P(i) 



+ 



and 



P(2) 



+ 



(D.9) 



The above exponents do not satisfy the quadratic (zero-mass-shell) Kasner relation, since 



X^G° M ) 2 — yZ~2iiP^J = — y • m addition, = !> which means that each metric's 

volume Vi expands as a function of its respective proper time as v oc t^, with i = 1,2. 
Note that the latter volume expansion has an accelerating behaviour and is exactly the 
same as in the p = case. The evolution is still anisotropic but slightly less than the 
Kasner case: 

2 



A x = A 2 



(D.10) 
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Let us now check, as we did for the p = case, if the obtained solution is stable against 
perturbations of 5. We have to recalculate the quantities A and B of Eq.(|8"9"D and their 
limits for r ^> 1. For r C both signs of p have: 



A = 



B 



mH 2 

< 

12 



,a/2 



(D.ll) 



and hence we have 5 = and absolutely stable motion. On the other hand, for the case 
t — > 0, since the evolution is exactly the same as for the p = case, the solution is unstable 
at late times as described in the main text. This instability is again linked to a run away 
of 7 towards large values. 
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